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The magnetic behavior of antiferromagnetic spin 1/2 chains with site-centered impurities in a 
magnetic field is investigated. The effect of impurities is implemented by considering different 
situations of both diagonal and off-diagonal disorder. The resulting magnetization curves present a 
wide variety of plateaux, whose position strongly depends on the kind of disorder considered. The 
relevance of these results to experimental situations is also discussed. 

PACS numbers: 75.10.Jm, 75.10.Nr, 75.60.Ej. 
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I. INTRODUCTION 

The study of the magnetic properties of low dimen- 
sional antiferromagnets has received much renewed at- 
tention during the last years. One particular issue that 
captured both experimental and theoretical efforts is the 
appearance of magnetization plateaux in spin chains and 
ladders systems. In general, those plateaux appear for 
pure spin systems at rational magnetization values [Q-Q ■ 
Some experiments have indeed confirmed these theoreti- 
cal predictions in a few particular cases || . 

More recently, the effect of the presence of impuri- 
ties on such magnetic behavior was also investigated. 
From the theoretical point of view, the properties at zero 
magnetic field have largely been elucidated HQ , and re- 
cent analysis explored also the robustness of the ratio- 
nal plateaux for small amplitude of the disorder ||. A 
very interesting phenomenon recently discovered was also 
that, for bond-disorder with discrete probability distri- 
butions, plateaux at non-rational values of the magne- 
tization are present ||. Moreover, the position of such 
plateaux in the curve is related in a simple way to the 
concentration of impure bonds in the system. Since this 
kind of disorder is a good candidate for modelling con- 
crete experimental situations, this result opens new per- 
spectives in the interpretation of experimental curves. In 
order to make closer contact with experiments, one has 
to take into account all the effects that can be produced 
by non-magnetic impurities, and provide a wide range of 
possibilities for implementing the presence of such impu- 
rities in realistic models. 

In this paper we analyze a simple model where different 
kinds of diagonal and off-diagonal disorder are present. 
The techniques used are based on a decimation proce- 
dure |Tc| |, as well as standard quantum mechanics tools 
and numerical calculations for XX chains of sizes rang- 
ing from 50 000 to 100 000 sites. Despite the simplicity 
of our model, we find that a huge number of families of 
cases is present, each of them presenting its own hierar- 
chy of magnetization plateaux. Moreover, we show how 
a simple tuning of the microscopic parameters can drive 



our system from one family to another. We also argue 
that this simple example serves as a good laboratory for 
understanding and classifying the wide variety of cases 
one can encounter in experimental curves of antiferro- 
magnetic systems with impurities. 



II. THE MODEL 

The model we are going to study is a spin- 1/2 chain 
whose Hamiltonian is given by: 

H = Y, [Ji (SfSf +1 + S«S? +1 + AS*S* +1 ) + hi Sf] .(1) 

i 

The randomness is implemented by considering Ji and 
hi as random variables. Specifically, bond randomness 
will be obtained by assigning to all the variables hi the 
same value fixed at h, and by giving a probability dis- 
tribution to the set of variables {Ji}- The probability 
distribution can, of course, be chosen in such a way to 
give, on average, a periodicity of q to the bond variables. 
It is important to say that realistic distributions relevant 
for many possible experimental situations are discrete. 
In this sense, a binary distribution captures the essen- 
tial characteristics of the phenomena we want to describe 
here, the generalization to more complicated discrete dis- 
tributions being straightforward. In |J, a distribution of 
the form: 



P(Ji) = p 5{J, - J 1 ) + (1 -p) 5{J, - J - 7i J) 



(2) 



was considered, where 7$ = 7, (—7) if i = qn, (i 7^ qn). 
In order to match the behavior of a disorder originating 
from a site impurity, a most appropriate distribution can 
be given by the following algorithm: distribute first reg- 
ular values for the bonds Ji with the desired periodicity 
q, and parameter 7, then, to each site of the chain assign 
a probability 1 — p to be an "original" site, and p to be an 
"impure" site. Finally, for each site i which turned out to 
be impure, change the values of the surrounding bonds 
Ji and Ji-x to J'. We are, as in concerned with three 
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different values for the bond strength, J±-fJ and J', but 
now with a correlated probability distribution. 

For site-centered disorder, we are going to consider two 
different cases. The first is just obtained by writing hi 
as h + h[ where h[ is a random variable taking values of 
with probability p and h' with probability 1 — p. This 
case can be considered as academic, since it violates the 
symmetry h — > —h, but it will nevertheless provide useful 
insight for more realistic disorders. The second and more 
realistic case is the h — > —h symmetry preserving case 
where hi = h(l + a^), with values for cij being binary 
distributed among and a with probability p and 1 — p, 
and a > — 1 (the limiting case of a = — 1 corresponding 
to an impurity which does not couple to the magnetic 
field). 

In the case of A = 0, the model can be mapped, via the 
well known Jordan- Wigner transformation, to a problem 
of free fermions, whose first quantized Schrodinger-like 
equation is: 



Ji—X — 1 ~\~ Ji hi Ci 



(3) 



The magnetization is simply related to the number of 
states occupied and the z component of the susceptibility 
is just proportional to the density of states. For a given 
energy, the number of states can be simply obtained by 
the node counting method JllJ , related to the number of 
positive "self-energy" variables: 

A; = Ci^i Ji-i/Ci. (4) 

The recursive formula for these variables is: 



A i+1 = J? /(hi - A,) 



(•5) 



which, for the case of bond and h — > —h symmetry pre- 
serving site disorder can be written as: 



A+i = Wf/{h - Di 



(0) 



where D, = 



and W 2 



, — jy-. — rj*-. r . Note that 

l+ai I (l+cti)(l+a;+i) 

Di has the same sign as At, then the node counting can 
be done in the D variables. This remark will be very 
useful when discussing the presence of plateaux and the 
behavior of the magnetization curve close to h = 0. 



III. OFF-DIAGONAL DISORDER 

As mentioned above, this case corresponds to a gener- 
alization for bond-disorder considered in . In this case 
each time a site is considered as an impurity, it must be 
surrounded by bonds with lower values. The decimation 
procedure used in joj is again well adapted for identifying 
the location of magnetization plateaux. However, some 
care has to be taken in counting the number of bonds for 
each decimation step, particularly in the case of two im- 
purities sitting in neighboring sites. We refer to the ref- 
erence mentioned above for the details of this reasoning 
and present here the result for a dimerised chain which 



is the most relevant case for experimental situations, the 
generalization to generic periodicity of the lattice being 
straightforward. The principal plateau is located at: 



M = 2p-p 2 



(7) 



coming from the decimation of the bonds with the high- 
est value. The second step of the decimation gives a 
secondary plateau at: 



M = 2p - 2p 2 + 2p A 



P 



(8) 



These plateaux are well observed in the numerical curves 
of figure (0), where we show the result for an XX chain. 

The behavior of the magnetization for small magnetic 
fields can also be obtained as a generalization of the nor- 
mal off-diagonal disorder. For example, for the XX case, 
a mapping to a random walk in the self-energy variables 
shows, as for a standard dimerised chain with random 
bonds |]l2|,0 , a power law behavior for the magnetic sus- 
ceptibility for even periodicity of the lattice: 



Xz^H 



A-l 



(9) 



where the exponent can be easily obtained from the mean 
and variance of In (J) (see for example ||). Using the 
same method, a logarithmic behavior for the susceptibil- 
ity is obtained for odd periodicity of the lattice: 



1 



Xz oc 



H [ln(ff 2 )f 



(10) 



IV. DIAGONAL DISORDER 

We first concentrate on the naive diagonal disorder one 
can introduce in a dimerised chain, by supposing that 
a supplementary magnetic field h' at each impurity is 
present. While for low h' no noticeable changes at the 
curve occur, for strong enough values of h', we see in fig- 
ure (^|) the appearance of new plateaux. A simple way 
to understand such plateaux is by considering the strong 
coupling case 7 ~ 1. The order zero in powers of (1 — 7) is 
just given by a combination of dimers which can contain 
two, one or zero impurities with probability (1 — p) 2 /2, 
p(l —p) and p 2 /2 respectively. It is then easy to draw the 
magnetization curve for each case (a three step stair like 
curve) and then superpose the curves with the appropri- 
ate weight. The result is compared to the numerical data 
in figure (J3J) . The positions of the plateaux predicted at 
this order are located at: 



M 



-1+p 2 , -(1-p) 2 , 0, p 2 , 2p-p 2 



fill 



Note in passing the lack of M <-> — M symmetry in 
the magnetization curve due to the very nature of the 
disorder. Comparing with the numerical result, one sees 
that some plateaux are softened due to the presence of 
non-zero coupling which was neglected in our procedure. 
The next step is as usual to turn on a non-zero value 
for 7 and use standard quantum mechanics perturbation 
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theory (see for example jHj for a similar treatment) . The 
result is a smoothing of the jumps between the plateaux 
and a correction to their width which can be calculated 
in powers of (1 — 7). Since we are concerned with more 
realistic kinds of disorder these computations are beyond 
the scope of this work. We just mention the fact that this 
case leaves the M = plateau present, a result which can 
be also obtained for the XX case by the random walk 
argument. We refer to Jll| for the details of this map- 
ping, and just mention here that it is easy to study the 
cycles of H — self energy variables for low energies. One 
can indeed show that, in the equivalent first quantized 
Schrodinger picture, and for small enough but non-zero 
values of the energy, these cycles are infinite indicating 
that the number of states remains constant in a finite in- 
terval of energy, which translate in the magnetic language 
into a zero magnetization plateau. 

The second case is obtained by assuming that each 
impurity couples to the magnetic field with a factor of 
(1 + a). It gives a symmetric magnetization curve by 
the change h — > — h. It is sufficient to study only the 
case a < 0. Indeed, the case of positive a is obtained by 
a simple re-scaling of the magnetic field and the change 
p <-> 1 — p. In the case of XX chains, we can make 
use of the change of variables defined in ([}]). We have 
then the equivalence between the model we are consider- 
ing and a system with constant magnetic field, but with 
a wide hierarchy of values for the bonds. The decima- 
tion procedure can then be used to predict the presence 
of plateaux. We have however to identify different cases 
corresponding to each hierarchic structure of the bonds, 
and these are given by the values of a and 7 we are con- 
sidering. As already stated before, the decimation proce- 
dure is implemented by eliminating the strongest bonds, 
and their relative weight in the system gives the position 
of the corresponding plateau in magnetization. We omit 
the details of the combinatorics and present the results 
for the position of the principal plateaux for each case in 
table (@). Figure (|) shows also the magnetization curve 
obtained for one of these cases, showing a fair agreement 
with the decimation predictions. 

Let us now analyze the behavior of the system at low 
magnetic fields. By using again the mapping to a ran- 
dom walk of the generalized Riccati variables, one could 
naively say that equation (Q) gives rise to a behavior 
similar to the one obtained for an off-diagonal disorder. 
An important difference is however that the noise in this 
case is cancelled step by step, giving rise to a regular 
and non-random walk at large times (this can be easily 
seen by inspecting the form of the variables Wi). Then, 
at low magnetic fields, the behavior of our systems is 
asymptotically equivalent to the one of a pure chain, and 
for an even periodicity of the bonds shares in particular 
the presence of a plateau at zero magnetization. Strictly 
speaking, these results are valid for XX chains only, since 
we have used the change of variables (g). We conjecture 
however that the qualitative behavior for the curve, and 
in particular, the position of the magnetization plateaux 



remains valid for generic XX Z chains. 

V. MIXED DISORDER 

This last case can be considered as the combination of 
both kinds of disorder studied before. It is important to 
notice that the equivalence between positive and negative 
a is not valid anymore due to the extra presence of bond 
disorder. For the XX case, it is however easy to see by 
means of (^|) and the subsequent decimation procedure, 
that the case a > has the same plateaux as for the off- 
diagonal disorder discussed above. This is still the case 
for negative but small values of a, since the hierarchy 
of coupling in the effective system is unchanged. The 
situation is however radically changed when a goes below 
the critical value: 



where the hierarchy of values for the bonds is completely 
changed. Figure (5) shows that crossing the critical value 
of a indeed changes completely the nature of the mag- 
netization curve and the position of the plateaux in the 
curve. In that case, the strongest bonds have the values 
J'/(l + a) and can be found in arbitrary long chains of 
impurities. This fact invalidates the standard decima- 
tion procedure for locating the position of the plateaux. 
However, one can proceed by noticing that, for values of 
a close enough to —1, the inter- impurities bonds in the 
effective model are much larger than the others. One can 
then use a kind of strong-coupling treatment in which the 
zeroth order system is obtained by equating the value of 
the remaining bonds to zero. We are then left with a 
collection of decoupled spin chains with arbitrary length. 
For an array containing n spins, the eigenvalues of the 
Hamiltonian are simply given by the roots of second or- 
der Chebyshev polynomials 

p7T 

E p = J ef f cos( — — ) ; p = l...n (13) 

71+1 

and J e ff is the value of the coupling between the spins. 
The magnetization curve for such a finite chain can also 
be easily obtained. To draw the magnetization curve for 
the total system one just has to superpose the curves for 
all possible values of n with the corresponding probability 
of appearance (1— p) 2 p n . The result is shown in figure @ 
where we compare the numerical result with the strong 
coupling computation performed up to 11th order in p. 
The appearance of a hierarchy of plateaux close to sat- 
uration predicted by the theory is clearly observable in 
the numerical data. 

The plateau at low values of the magnetization and the 
jump in the magnetization at zero field is an artifact of 
the zeroth order approximation. This can be cured now 
by turning on the remaining coupling to non-zero val- 
ues and treating them in perturbation theory. On what 
concerns the jump at M — 0, there is however a much 
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simpler way of studying the shape of the curve close to 
Al = 0. This is achieved again by the mapping to the ran- 
dom walk problem. One recovers in this way the power 
law or logarithmic behavior as in (p|)-([l0|) depending on 
the periodicity of the bonds. 

VI. CONCLUSION 

To summarize, we have shown that the magnetization 
curve of spin chains with different kinds of impurities ex- 
hibits a rich variety of plateaux. While the phenomenon 
of non-rational plateaux was already observed in a dis- 
ordered case ||, we have shown in this paper with sim- 
ple examples that the position of such plateaux strongly 
depends on the particular values of the microscopic pa- 
rameters. While one can speak of a sort of universality 
for the position of the plateaux, depending only on the 
concentration of impurities and periodicity of the aver- 
aged bonds, we have clearly seen that different families 
of disorder produce qualitatively different magnetization 
curves. Moreover, within the same kind of disorder (the 
mixed disorder case in particular) a simple tuning of the 
microscopic parameters can switch the system from one 
"universality class" to another. This scenario teaches us 
that the characterization of spin chains with impurities is 
better understood in terms of families of disorder, at least 
to understand the behavior of the magnetization curve. 
While some results in this paper have been obtained for 
the XX case only, it is reasonable to think that such a 
characterization in terms of families of disordered chains 
remains valid for generic XX Z chains. Of particular in- 
terest is the mixed case, for which we believe the case 
a — y — 1 should represent the magnetization curve for 
materials such as CuGe03 doped with Si |Q. We trust 
this work provides the necessary tools for predicting the 
shape of the magnetization curve in future experimen- 
tal situations, where higher periodicities than 2 are also 
conceivable 

We acknowledge B. Dujardin for fruitful discussions 
and D. C. Cabra and A. Honecker for constructive com- 
ments and a careful reading of the manuscript. 
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FIG. 2. Magnetization curves for the first kind of diago- FIG. 4. Magnetization curves for the (h <-> -h) symmetric 

nal disorder mentioned in the text (see eq. ( |ll| ) taken with disorder with J = 2, a — —0.7, 7 = 0.7. 
p = 0.3, 7 = 0.5, J = 2. 
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FIG. 6. Strong-coupling and numerical result for the mixed 
disorder case with J — 2, J' — 0.8, a — —0.95, 7 = 0.5, 
p = 0.3. 
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TABLE I. Magnetization at the principal plateaux for the 
(h <-> — h) symmetric diagonal disorder (see text). 
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